Quantum pumping in a ballistic graphene bilayer 
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We investigate quantum pumping of massless Dirac fermions in an ideal (impurity free) double 
layer of graphene. The pumped current is generated by adiabatic variation of two gate voltages in 
the contact regions to a weakly doped double graphene sheet. At the Dirac point and for a wide 
bilayer with width W ^> length L, we find that the pumped current scales linearly with the interlayer 
coupling length l± for L/l± <C 1, is maximal for L/l± ~ 1, and crosses over to a \n(L/l±) /(L/l±)- 
dependence for L/l± S> 1. We compare our results with the behavior of the conductance in the 
same system and discuss their experimental feasibility. 

PACS numbers: 73.23.-b, 72.80.Vp, 73.23.Ad 



I. INTRODUCTION 

Quantum pumping of charge refers to the generation 
of a dc electrical current in the absence of an applied bias 
voltage by periodic (ac) modulation of two or more sys- 
tem parameters, for example the shape of the confining 
potential or a magnetic fieldi. The idea of adiabatically 
generating a flow of particles in a moving periodic po- 
tential is due to Thouless 2 , and has been followed by 
many theoretical and a few experimental investigations 
of pumping in mesoscopic systems. In 1998 Brouwcr 3 , 
building on earlier results by Biittikcr et at— and a pro- 
posal by Spivak et alfi, developed a description of quan- 
tum pumping through open mesoscopic systems in terms 
of the scattering matrix of the system. This paved the 
way for investigating the effects of quantum interference 
on quantum pumping and has led to theoretical investi- 
gations of many different aspects of pumped currents in 
open nanodevices, such as the relation of quantum pump- 
ing to geometric (Berry) phases^, the effect of Andreev 
reflection on quantum pumping in nanostructures that 
contain superconducting partsi, the effect of electron- 
electron interactions^, and the generation of adiabatically 
pumped spin currents^. 

Most of these investigations were carried out for 
semiconductor structures such as quantum dots and 
nanowires. In addition, also proposals for adia- 
batic pumping in carbon nanotubes 10 , and recently 
graphen e 11 ' 12 have been put forward. Since its exper- 
imental discovery in 2004^2, graphene has been found 
to exhibit electronic transport properties that are quite 
different from those in other nanoclcctronic structures 
due to the nature of its charge carriers (massless Dirac 
fermions described by a relativistic wave equation) 1 —. An 
example of this is the importance of evanescent modes 
for transport close to the Dirac point: in a sample of 
undoped graphene, which does not have any free elec- 
trons, contacted by doped electrodes the conductance 
close to the Dirac point is dominated by the contribu- 
tion of evanescent modes which transmit electrons in- 
jected from one end of the sample to the other en d 15 ' 16 . 
This is also true for quantum pumping in a monolayer of 
graphene, where the pumped current is induced by two 



oscillating gate voltages^. 

In this paper we investigate quantum pumping of Dirac 
fermions in a graphene bilayer. Compared to a carbon 
monolayer, the bilayer has an additional energy scale, 
namely the interlayer coupling strength t±. The corre- 
sponding lengthscale is an order of magnitude larger 
than the interatomic distance S^. Our aim is to investi- 
gate the dependence of the adiabatically pumped current 
on the interlayer coupling t±. For the conductance in a 
bilayer with heavily doped contact regions and width W 
3> length L it has recently been found that at the Dirac 
point the bilayer transmits as two monolayers in paral- 
lel and the conductance is independent of t±— . In con- 
trast, for the pumped current I p we find that I p depends 
linearly on L/l±_ for small interlayer coupling strength 
t± = 1/Zj_ <C 1/L, exhibits a maximum around <j_ ~ 1/L 
and scales as \n(L/l±)/(L/l±) for large interlayer cou- 
pling ij_ = 1/L. For typical experimental pa- 
rameters the pumped current is of order 10-100 pA. 

The paper is organized as follows. In Sec.[TT]we present 
the bilayer model and an introduction to pumped cur- 
rents. Sec. IIIII contains the calculation of the pumped 
current, followed by results and comparison to the con- 
ductance in Sec. IIVI We conclude by making a connec- 
tion to experiments in Sec. |Vl 



II. BILAYER MODEL AND PUMPED 
CURRENT 

We consider the geometry that is schematically de- 
picted in Fig. [1] A sheet of ballistic graphene in the 
(x, y)-plane contains a weakly doped strip of length L 
and width W which is contacted by two more heavily 
doped electrodes at x = and x — L. The doping in 
these contacts is controlled by gate voltages, which in- 
duce a potential profile of the form 

U\(t) for x < or x > L in the upper layer 
U (x) = < U%{t) for x < or x > L in the lower layer 
for < x < L in both layers 

(1) 
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FIG. 1: Schematic picture of the graphene bilayer. Top panel: 
Two stacked honeycomb lattices of carbon atoms in a strip 
of width W between metallic contacts (blue and red regions). 
Bottom panel: Variation of the electrostatic potential across 
the two layers. 



We assume the potential step to be abrupt, which is 
justified close to the Dirac point where the Fermi wave- 
length Xp > L and any smoothing of the step over a 
distance small compared to L becomes negligible. In 
addition, we consider a short and wide geometry (L <C 
W), for which boundary conditions in the transverse y- 
direction become irrelevant. 

The bilayer pump is operated by periodic variations of 
the carrier density in the leads by varying the potentials 
Ux(t) = Ui + 5Ui cos(wi) and U 2 (t) = U 2 +SU 2 cos(ut+(j)) 
such that U\ (t) = U 2 (t) on average. In the linear response 
regime where SUi « f/j (i = 1,2), the pumped current 
I p into the left lead is given by the scattering matrix 
expression 3 



cue svn.<j>8U\8U-> v— * T f d^aB dS a R 

^ = s 2^ 2^ Im 



2 it 



V dUi dU 2 



Here the index a sums over all modes in the left con- 
tact region and f3 sums over all modes in both the left 
and right contact regions. The pumped current into the 
right lead is then given by Ir = — 7j, = — I p . S denotes 
the Landauer-Biittiker scattering matrix whose elements 
S a /3,nm describe scattering from mode m in lead j3 to 
mode n in lead a. 

In the presence of a potential U\ (U 2 ) in the upper 
(lower) layer, the low-energy excitations of the graphene 
bilayer close to a Dirac point are descibed by the 4x4 
Hamiltonian, H = 



H= 



t 







( Ui v(p x + ip y ) 
v(p x - ip y ) Ui 

i_L U 2 v(p x - ip y ) 

v(p x + ip y ) U 2 



,(3) 



where p = —ihd/dr is the momentum operator. The 
Hamiltonian (J3|) acts on the four-component wavefunc- 
tion (ipAii ipBt, ipB 2 i i J A 2 )> with A\ labeling the ampli- 
tude on the A-sublattice of the first (upper) layer, and 
similarly for Bi, A 2 and B 2 . We only take nearest- 
neighbour coupling from A to B sites within the same 



layer or between the two different layers into account^. 
The four eigenenergies of the Hamiltonian (J3J are given 
by 



£1,2 



U!+U 2 



U1 + U2 1 rriTK 

£3,4 = ^ 



(4) 



where f±(k) 



4fc 2 



2^/t 4 ± +4k 2 [t 2 ± 



2t\ + {Ui - U 2 ) 2 ± 

(C/i - U 2 ) 2 ], and k = (kl + fc2)i/2 de _ 
notes the total momentum. From now onwards, we ab- 
sorb a factor (hv) -1 in e, Ui, U 2 and t±, which are all 
given in units of momentum. 

Using scattering matrix theory we calculate in the next 
two sections the total pumped current 1 - I p for U\ = U 2 = 
U at the Dirac point e = and derive analytic expressions 
for the limit of heavily doped contacts U — > — 00 (these 
were also considered in RefJ^). 



III. CALCULATIONS 

The scattering matrix S and subsequently its deriva- 
tives with respect to U x and U 2 can be found by matching 
eigenstates of the Hamiltonian ([3]) at the interfaces x = 
and x = L. For given e and transverse momentum k y the 
eigenstates of are characterized by two longitudinal 
momenta k x ± 



V? 



V 2 



± =:d{V? - V 2 ) 2 + U\ViV 2 - k 2 



(5) 

with Vj = e — Uj (j — 1,2). Associated with each real 
wavevector k x+ are two propagating modes <f>^Ax,y) 
(right-going) and <p^ + (x,y) (left-going). Similarly, an- 
other two propagating modes (j)f_(x,y) and 4>^_(x,y) 
correspond to each real wavevector k x -. Defining k± 
k . . + iky, the left- and right-going eigenstates are given 
by 



(x,y) = N± 




-ik^-i-x+ikyy 



where 
X j± 



kl± 



k 2 + t ± V 3 



1,2 



N± = [2Wk x ±(V x Xl ± + V 2 Xl ± )] 



-1/2 



(6a) 



(6b) 



(7) 
(8) 



The eigenstates ^ are normalized by N± such that each 
state carries unit current 



I 

ev 



dy<j> 



a x 
a x 



= 1. 



(9) 



oo Eqns. ([5]) and (J6]) reduce to the integrating over k y then yields the pumped current: 



For U x = U 2 - 
results of Rcf. 17 . 

From now onwards we assume to be at the Dirac point 
e = 0. In the undoped graphene region < x < L 
we then find from Eq. ([5]) that k x ± = ±ik yi which cor- 
responds to evanescent modes. The left-incident eigen- 
states of the Hamiltonian ([3]) can then be written as 



</>? t± {x,y)+ rXt <p%+ {x,y) + r L _ L ± <%_{x, y) 



[(ci±Xi +c 2 ± x 2 )e fc ^+ 
(c 3 ± X3 + c 4 ± Xi) e ~ ky 



for x < 



SI e ik y y 

for < x < L 
>*-(x-L,y) 
for x > L 

(10) 

denote the 
mode 

incident from lead /3 to a — mode in lead a. In the mid- 
dle region the eigenvectors Xi"X4 have been constructed 
such that they are linearly independent at the Dirac point 
e = 0. This yields 



Here 



a/3 



and t 



a/3 



with 



_ _+ ^ rr a,0 e {L, R} 

reflection- and transmission coefficients from a 



Xi 



Xa 




(11) 



The reflection and transmission coefficients are calculated 
by matching the eigenstates (JTDJ) at the two interfaces 
x = and x = L, see Appendix [X] For a short and 
wide geometry with L W the boundary conditions in 
the y-direction become irrelevant. Taking infinite mass 
boundary conditions, such that k y is quantized as k y 
: n. + l/2)n/W, n = 0, 1, 2, . . . the pumped current fl5J 
becomes a sum of eight terms: 



with 



lu e W sin 8 U\ 8 U 2 
2^2 



dk y J2 u ^'( k y) (12) 



LL 



d& at 



LR 

go' 



\ dVx dV 2 dVx dV 2 



IV. RESULTS 



In Appendix [5] we describe details of the calculation 
of the derivatives of r^, and t^, with respect to U\ 
and U 2 and consider the limit Ux -+ U 2 = U with [/ — > 
This is equivalent to Vx — >• V 2 = V 7 with V — > 



-00. 



-00. Substitution of these derivatives into Eq. (TT2"j) and 



/p = 



weWsin05J7i(H72 4A 



dkaj 



2tt 2 c/2 y 

(A 2 +6fc a Lsinh(2fc y i)-8cosh 2 (fc a i)) sinh(2fc a L) 



(A 2 + Acosh 2 (k y L)Y 



2uje W 



A sine/) 



8U X 8U 2 
U 2 



VL-s-oo 



dx 



(A 2 + 6xsinh(2a;) - 8 cosh 2 (x)) sinh(2x) 

(A 2 +4 cosh 2 (a;)) 3 
uje W . 8U X 8U 2 1 

L Sm( ^ C/ 2 4A 2 (A 2 +4) 2 



A 5 - 20A 3 + 3A(A 2 + 2) (A 2 + 4) ln(A 2 + 4) 



(14a) 



(14b) 



+12^2+4^2(2 + y + ^VA 2 +4) 
A 2 A 



-I/ J (2 + y-- % /A 2 + 4 



(14c) 



Here we have defined the dimensionless coupling length 
A = t±L and Lz 2 is the dilogarithm function defined as 
Li 2 = f*Qnt)/(l - t)dt. Eq. (HU is the main result 
of this paper. The upper integration limit in Eq. (|14a[) 
originates from the requirement that the wavefunctions 
in the left and right leads should correspond to travel- 
ing waves, and hence that k x ± should be real (k x - gives 
the more stringent condition k y < y/V(V — t) « V for 
V — > 00). We see that both for L = and for L -+ 00 
the pumped current (|14c|) reduces to zero. For L = 0, 
i.e. in the absence of the middle region, there is no possi- 
bility for evanescent waves to interfere in this region: all 
incoming waves from one lead are fully transmitted into 
the other lead without scattering to the other layer and 
hence do not contribute to the pumped current (which 
is composed of waves that scatter at least once from one 
layer to the other). For L — > 00, the evanescent modes 
in the weakly-doped middle layer do not reach the other 
contact and the pumped current becomes zero. 

Fig. H shows E CT , CT '=± IW(fc y ) [Eq. (lATTT) in units 
of 4A?7~ 2 ], which is essentially the mode- dependent 
pumped current, as a function of k y . We see that 
X^ct<t'=± ^tra-'(k y = 0) = 0, i.e. waves with transverse 
(13) momentum k y = that are incident perpendicular to the 
interface do not contribute to the pumped current, which 
is a manifestation of the Klein parado x 14 ' 15 . We also see 
that negative mode-contributions only occur for small 
values of A. From Eq. (|A1 1|) it follows that the largest 
contribution to „i—± n acr r comes from the transverse 
modes k y L ~ 0.1 A for A ~ 50, similarly as for the 
transmission (conductance) in the same system 1 !. Fig. [3] 
shows a plot of the total pumped current I p as a func- 
tion of A, for a fixed value of W/L = 100. For n ~ 200 
modes the sum converges to the continuum result. The 
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-0.06 




by dI p /dX = and yields 



3.E 



/ P (A n 



nci we If . .SUM 

= - 51 2^T sm ^-^' 

(17) 

Thus I p is maximal if ij_ is of the same order as 1/L. 
The pumped current thus strongly depends on t± and 
reduces to zero for t±_ — > 0. This is in sharp contrast with 
the behavior of the conductance at the Dirac point in a 
graphene bilayer, which equals the conductance across 
two monolayers and is independent of t±— . 



FIG. 2: The mode-dependent pumped current Eq. (|A11|) as a 
function of k y L. The curves correspond to A = 0.01 (purple), 
A = 2 (red), and A = 5 (blue). 



CONCLUSION 




FIG. 3: The total pumped current I p as a function of A, for 
W/L = 100. The blue, purple and red lines represent the 
exact sum over k y of the integrand in Eq. (|14a[) for n=100, 
150 and 200 modes. The black dotted line is the integrated 
result Eq. (flic)) . 



number of modes needed for convergence becomes larger 
for W/L larger, while for W/L = 20 only n ~ 20 modes 
are needed. As expected, I p (t± = 0) = 0, since in the 
absence of interlayer coupling the bilayer reduces to two 
uncoupled monolayers for which the pumping parame- 
ters Vi(t) and V-zit) become uncoupled and no pumping 
occurs. From Eq. (|14cl) we obtain that for small inter- 
layer coupling strength < 1/L the pumped current I p 
scales as 



lim In 



41n2 



64tt s 



1 W . 

— toe — sin i 



SU 1 SU 2 

u 2 



(15) 



hence I p depends linearly on t± for t± = l/l± -C 1/L. On 
the other hand, in the limit of large interlayer coupling 
tj_ = I fix > 1/L wc find 



In summary, we have investigated the adiabatically 
pumped current through a wide graphene bilayer con- 
sisting of a central weakly doped graphene sheet coupled 
to two heavily doped contact regions. At the Dirac point, 
the pumped current is carried by evanescent waves in the 
central region and exhibits a cross-over from linear (for 
t±L <C 1) to logarithmic (for txL ;§> Independence as 
a function of increasing interlayer coupling strength tx, 
with a maximum around ij_L ~ 3.88. This scaling be- 
havior with tx is markedly different from the behavior 
of the conductance G in the same system, which is inde- 
pendent of tx and equal to the conductance across two 
monolayers in parallel. In practice, this different behav- 
ior of Ip and G as a function of tx and L could be used 
to distinguish between the conductance and the pumped 
current. 

We can estimate the magnitude of the pumped cur- 
rent using typical experimental parameter a ' 20 tx = 
{0AeV)/(hv F ) w 0.6 10 9 m" 1 , U « 0.1 V, SU « 10 
mV, and u « 1 GHz. In order to be able to measure 
a substantial pumped current, one thus needs sample 
sizes L ~ 4/ij_ ~ 6 nm, which is smaller than currently 
available samples of order /jm^. However, with steady 
progress towards smaller sample sizes, sample lengths of 
order 10-100 nm are expected to become available in the 
future. We then obtain from Eq. (|14c[) : 



Ip,max ^ 5 ' 10 



(18) 



which, for W/L ~ 100 — 1000 is well within experimental 
reach. Observation of this current would be a striking 
demonstration of quantum pumping produced by rela- 
tivistic quantum mechanics. 



lim /„ 



3 W 
—r oje — sim 

47T 2 L 



SUM InA 
~lP A 



(16) 
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Appendix A: Calculation of the scattering matrix 
elements and their derivatives 



In this appendix we calculate the reflection and trans- 
mission coefficients r^, and t°f G , (a, a' 6 {+, — }; a, € 
L, R) from Eq. ([TO]) and their derivatives with respect to 
Ui and U 2 . These derivatives are then used to calculate 
the pumped current (fT2|) . 

Matching the eigenfunctions (fTO)) at the interfaces x = 
and x — L results in the 8 equations: 



and 



, LL 



LL 



a 

,,LL 



a +i (subindex + O subindex 
Pi— (subindex + O subindex 



= 7_i__ (subindex + O subindex ■ 

= (subindex + -f-> subindex 

= (subindex + <-> subindex 

7"~ = 7++ (subindex + O subindex ■ 



r + ±JV 4 



/ 


^2+^1 \ 




( x 2 _v x \ 




-x 2+ k+ 


+ r- ± N- 


-X 2 -k- 




x 1+ v 2 


X!-V 2 


\ 


-X 1+ kl J 




V -x^k*_ j 




N± 



I X 2 ±Vi 
X 2± k* ± 

x 1± v 2 

\ X 1± k± 



( X 2+ Vi 
X 2+ k* + 

x 1+ v 2 
V x 1+ k+ 



+ t_±JV_ 




/ c 3 ±z~ l 

(ci± - itj_Lc 2 ±)z 

V (c4± - i<j_ic 3 ±)z _] 



(Al) 



(A2) 



with z = exp(fcj,L), Vj — —Uj at the Dirac point 
(j = 1,2) and Xj± and A± given by Eqns. and ©. 
Eliminating ci — C4 from Eqns. ()A1[) and (|A2|1 yields af- 
ter some straightforward but lengthy algebra for the four 
reflection coefficients: 



Also 

5 = 2N+N- \—{k x+ k x - — ky)(A 2 + B 2 ) 

+ (k 2 x+ + k 2 x _ - 2k 2 )AB] (A6a) 
e = 2iN+N_ [-(k x+ + k x ^)k y (A 2 +B 2 ) 

+ 2 (k x+ + k x -)k y AB] (A6b) 
T) = N+N- [-\ 2 {A- B) 2 VxV 2 

-2iX(k x+ - k x .){A - B)(X 1+ X^V 2 - X 2+ X 2 -Vx) 

-2{k x+ k x - + k 2 y ){A-B) 2 -2{k x+ -k x J) 2 AB] (A6c) 

The transmission coefficients in Eq. (|A2|) are given by 



[rti cosh(fcj,L) — v 



rl , 



\Yl(kyL) 



(A7a) 



t RL _ M++ vA^uyvir' 

++ _ <5 cosh(2fc y L) + e sinh(2fcj / £) + 77 

t RL = V+i cosh(k y L) - v!*i smh(k y L) 



RL 

-+ 
RL 



5 cosh(2k y L) + esmh(2k y L) + r\ 
t+i (subindex + <H- subindex -) (A7c) 
fc5+ (subindex + O subindex -). (A7d) 



with 



,,LL 



.ll _ cosh(2k y L) + gffi smh(2k y L) + 7 —. 
(5 cosh(2fcj,L) + e smh(2k y L) + rj 



(A3) 



with a, a' G {+,—}. Using k± = k x ± + ik y , the di- 
mensionless coupling length A = tj_L, A = X\ + X 2 _ and 
B = Xi-X 2+ , the coefficients in Eq. (|A3|) are given by: 



7^ 



:) +- 
„LL 



-- 2N + N_ [ik y {k+A 2 - k* + B 2 ) + 

= 2JV+ N-k x - [k+A 2 - k\B 2 - 2ik y AB] 
-- N + N_ [\ 2 {A- B) 2 ViV 2 - 
2iX(A-B)k x -(X 1+ X 1 _V 2 -X 2+ X 2 _V 1 )} - 
-- 2 N 2 h x _X 1 _X 2 _ [{k x+ - h x _ + 2ik y )A 

+ (k x+ - k x - - 2ik v )B] 
-- 2N 2 _k x -X 1 -X 2 -{k, x+ + k x -){A-B) 
-- ~2i\N 2 _ (A- B)k x _(X 2 _V 2 - X 2 _Vx) - 



LL 



(A4) 



a 



LL 



A*++ = -2 N+ A_ (A - B) k x+ [2k x -(A - B) 
+ i\(X 1+ X 1 _V 2 - A-2+AVFx)] 
= Ai N+ A_ (A - B) 2 k x+ k y 
fi+i = —2i A N 2 (A - B) k x - (Xf_ V 2 - X\_ Vi ) 
v rl = 2N_(A-B) k x - [i\(X 2 _V 2 + X|_Fi) 

-2X^X 2 ^{k x+ -k x ^)} (A8) 

In two limiting cases the reflection- and transmission co- 
efficients (|A3|) and (|A7I) reduce to simple forms : 

1) For L — > 0, the +- and — modes decouple and we 
obtain from Eq. (| A7|) 



t++ 1 and -> 0. 



(A9) 



2) For fj_ — ^ 0, i.e. in the absence of interlayer cou- 
pling, the reflection- and transmission coefficients (|A3[) 



6 



and (| A7|) reduce to the monolayer expression: 



,15,16 



+ + 

|yi|,|v 2 |-*-cx3 



— iky) sinh(fcj,L) 
fc x _l_ cosh(/c a L) + iky sinh(fcj,L 



J.RL 



|Vl|,|V2|- 

-4 



tanh(fcj,L) 

1 1 — fc. 



(AlOa) 



k x + cosh(k y L) + ifcy sinh(fcj,L) 

' (AlOb) 



cosh(fc y L) 



(AlOc) 



The reflection- and transmission coefficients for right- 
incident Dirac fermions rfffi and t^, ((J, a' G {+, — }) 
are given by Eqns. (|A3|) and (| A7|) by interchanging the 
layer subindices f and 2. So r++ = 7++ (subindex 1 o 
subindex 2) etc. 

In order to calculate the pumped current (|12l) . we 
need the derivatives of r^, and t^p with respect to Z7i 
and t/2- These can be calculated by splitting the coef- 
ficients orri etc. in Eq. (|A3|) into real and imaginary 
parts and differentiating each of these with respect to 
V\ and V2. The latter does not make any difference for 
the current, since only products of derivatives with re- 
spect to U\ and U2 enter Eq. (fT2|) which are the same 
as those of the corresponding derivatives with respect 
to V\ and V2 (since Vj — —Uj for e = 0, j = 1,2). 



The resulting expressions for the derivatives are rather 
lengthy and therefore not given here. For equal bias 
voltages V\ — V2 = V (as required for a true pumping 
process) and in the limit of V — > 00, however, an ana- 
lytic expression for the pumped current I p [Eq. |(5J) can 
be derived. This expression is obtained by first taking 
the limit V\ — > V 2 , for which X\- — > X 2 - = 2t±V and 
X 1+ -> -X 2+ = (V 1 -V 2 )(V+±) ->0 (see Eq. ©). We 
then expand all terms in the Eq. (|T3t in orders of V\ — V 2 
(or, equivalently, A 1+ ) and find that while the numer- 
ators contain terms of order (V\ — V2) 7 and higher, all 
terms in the denominator are proportional to (Vi — V2) 8 . 
In the limit of V — > 00, assuming V 3> t±,k y and 
Vt± ^> ky, the terms of order (V\ — V2) 7 in the numer- 
ator cancel and we are left with a finite pumped current 
which to highest order in V amounts to: 



4A 



^ ^-cra'(ky) — 



<7,<t'— ± 



(X 2 +6k y L s'mh(2k y L) -8 cosh 2 (k y L)) sinh(2fc y L) 



(A 2 +4cosh 2 (fc y i)) ; 



(All) 



Integration of (|A1I[) over the transverse modes k v then 
leads to the pumped current Eq. (|f 4[) in the main text. 
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